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ABSTRACT 


The forces acting on yawed smooth and rough circular 
cylinders in oscillating flow have been investigated for 
the purpose of determining the appropriate force-transfer 
coefficients and the applicability of the "independence 
principle." The results have shown that the flow about 
each cylinder is unique and the independence principle 
does not hold true. In addition, a comprehensive sensi- 
tivity analysis has been carried out to investigate the 
effects of undetected environmental disturbances on the 
force-transfer coefficients and the limits of probable 


errors have been delineated, 
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The study of time dependent flow about yawed cylinders is 
a topic which presently has great interest, both theoretically 
and practically. Of the great variety of time-dependent flows 
about yawed cylinders, the wave motion and Sinusoidally oscil- 
lating flow are of such importance as to warrant immediate 
investigation. For lack of adequate information concerning 
time-dependent flows, industry is often forced to adopt 
steady flow relations to time-dependent flow situations. 

This research proposes to investigate the forces acting upon 
yawed circular cylinders immersed in sinusoidally oscillat- 
ing uniform flow and to examine the validity or limitations 
of the present methods of analysis. 

im 1950, Morison, et al. {1] introduced an equation for 
calculating the in-line force due to unbroken surface waves 
acting on a vertical pile. For a cylinder of diameter D, 
the force per unit length is expressed as: 

2 DU 


mo eDee eur O-cop 4D Cos (1) 


bt? 


where U represents the incident flow velocity; Cas the drag 
coefficient; and Cat the inertia coefficient. The coeffi- 
clients Cou and Ca were considered to be time-invariant and 


constant along the length of the cylinder. Morison's paver 


i> 





was intended as a preliminary report with follow-on studies 
to be done on other structures in various wave actions. 
Morison, et al. did not consider the contributions of 
transverse forces and vortex shedding in the calculation 
of the in-line forces. 

Attempts have been made to extend this rather simple 
relation to far more complex situations such as combined wave 
and current flow, hydro-elastic oscillations, and yawed 
cylinders. Heideman, etal. [2] studied Morison's equation 
using ocean data and concluded that Morison's equation was 
Satisfactory for normal cylinders. The utilization of an 
equation unproven even under ideal flow situations to 
correlate ocean data so as to prove the validity of the 
equation is not a very meaningful exercise. Sarpkaya and 
Isaacson [3] reviewed in detail the numerous studies regard- 
ing the application of the Morison equation. 

It 1S obvious that there is a great need for an ideali- 
zation of the problem, or an experiment which is more 
manageable (e.g., a sinusoidally oscillating uniform flow). 
Only in this way can all the complex interactions be 
separately taken into account. 

Engineers faced with the problem of dealing with wave 
forces on yawed members, having no other recourse, drew 
upon previous work with steady flow. Hoerner [4] proposed 
the "independence principle", which stated that the normal 


pressure forces are independent of the tangential velocity 
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Bee SubCritical values of Rew where Re. is the Reynolds 
number based upon the component of the flow velocity normal 

to the cylinder. This principle allowed the neglecting of 

the tangential components. Bursnall and Loftin [5] found 

that the independence princivle does not apoly to the critical 
and transcritical flow regimes. Norton, et al. [6] found that 
the indevendence principle does apply to the postcritical as 
well as subcritical flow, but not to the critical and trans- 
critical regions in between. Thus, recent research has shown 
that the independence principle applies when the boundary 
layer is wholly laminar (Hoerner) or wholly turbulent (Norton), 
but its use in the critical and transcritical regions is 
uncertain. 

Tne designers of offshore structures were thus led to 
adopt the indevendence principle forthewave force calcula- 
tions and, in doing so, to generalize the Morison equation. 
This has been acceoted practice in industry for lack of a 
better relation. It can be asked what would be anticipated 
for oscillating flow or waves on the basis of what is known 
for steady flow. The instantaneous Reynolds number in such 
a flow will vary from SRG a to TR during a complete flow 
cycle. It could be postulated that the boundary layer would, 
at times, be fully laminar; at other times fully turbulent; 
and the rest of the time be in transition. In light of this, 
it is rather doubtful that the independence principle applies 


femall to oscillating or wave flow. 


ey 





The study of forces acting on yawed cylinders can be 
investigated either by oscillating a yawed cylinder ina 
tank, using small amplitude waves in a laboratory channel, 
@eerllating the flow about a fixed cylinder, or by using 
ocean data if available. Oscillating the cylinder has 
proven to be impractical because of difficulties in account- 
ing for the inertial force acting on the body, producing 
repeatable oscillations, controlling Spurious vibrations in 
the system, and in measuring the in-line and traverse forces 
Simultaneously. Waves have relatively more-complex flow 
kinematics due to the orbital motion of particles and the 
decay of wave amolitude with depth. Ocean data are not 
available for such a study. Oscillation of the flow past 
the cylinder has vroved to be the most practicable method. 

By using a sinusoidally oscillating uniform flow (which can 
be fairly easily produced) the kinematics of the flow are 
accurately known. Sarpkaya and Isaacson [3] discuss the 
Similarities and differences between such a flow and the 
Flows described by wave theory. 

The primary objective of this investigation was to study 
the forces exerted by a sinusoidally oscillating uniform flow 
on yawed circular cylinders to determine whether or not the 
independence principle is applicable. If the principle holds, 
the force transfer coefficients calculated by Fourier analysis 
using the normal velocity component should reduce identically 


to the normal cylinder case at corresponding values of K, Re, 


eS 





and k/D. If the independence principle does not apply, it 

is desired to determine what the coefficients are as functions 
of yaw angle, roughness ratio k/D, Re, and K. It would also 
be necessary to determine how well Morison's equation works 
with the new coefficients. 

With the foregoing objectives in mind a detailed investi- 
gation was undertaken using smooth and rough cylinders of 
nominal 6-inch, 45-inch, and 3-inch diameters at yaw angles 
of 45 degrees, 60 degrees, and 90 degrees (as measured between 
the ambient flow direction and the cylinder axis). Part of 
this effort has already been reported [7]. The present study 
extends the previous work, accomplishes the original objec- 
tives of the investigation, and takes a much closer look at 
the sensitivity of the modified Morison equation and of the 


Measured forces to random external disturbances encountered 


maenature. 
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II. EXPERIMENTAL APPARATUS AND CYLINDERS 


The experiments were carried out in the large oscillating- 
flow tunnel at the Naval Postgraduate School. The tunnel was 
basically unchanged from the configuration described in [8 
and 9]. Mounting pads for the strain gage housings were 
placed on opposite sides of the tunnel in such a manner so 
as to accommodate the yawed cylinders (Fig. 1). The mounting 
pad on one side of the tunnel was fitted with a slide so that 
the pad could be adjusted small distances to exactly match 
the cylinder length. The force transducers and housings were 
unchanged from those described in [9]. 

The cylinders also were similar to those previously 
reported, except that the ends were cut exactly at the 
desired angle of yaw. The cylinder ends were parallel to 
the tunnel walls with a 1/32-inch gap at each end. This gap 
was filled with a soft foamy material glued to the ends of 
the cylinder. Ball bearings were mounted in the ends of the 
cylinder with theouter bearing faces flush with the face of 
the cylinder. 

The cylinders mounted in the tunnel were supported by 
the pins in the force transducers. The pins were pervendicu- 
lar to the tunnel wall. Rubber washers were placed behind 
the pins so that the lateral motion of the cylinders (normal 


to the tunnel wall) was restrained, but without compressing 
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the cylinder. The cylinders were free to rotate, but a small 
clip was placed on one end of the cylinder to restrict the 
rotation to a very small amount so that the yawed end of the 
cylinder did not strike the tunnel wall as it rotated. Thus, 
as mounted, the cylinder was laterally constrained, but free 
to rotate within a small range. Careful calibration and 
extensive testing proved that this mounting system allowed 
accurate and repeatable recording of the normal and trans- 
verse forces acting on the cylinder. 

Calibration was first conducted in a vertical direction, 
Since yaw would not affect the total force in this direction. 
Weights were hung from the center of the horizontal cylinder, 
and calibration factors were obtained by converting the sum 
of the signals of two force transducers to pounds/mm and 
pounds/volt. The strain gages were found to be linear 
throughout the range of expected forces. The horizontal 
forces normal to the cylinder are related to the streamwise 
force by thesine of the yaw angle a (as measured between the 
cylinder axis and ambient flow). A system of supports and 
pulleys allowed applying horizontal force at the center of 
mvercylinder and normal to it (Fig. 1). The sum of the 
outputs of the two transducers was recorded on the strip 
chart recorder, establishing the necessary relationship 


between the normal in-line force* F and the total electrical 


* 
The force in the plane defined by the velocity vector and 
the cylinder axis. 


Zee 





output of the gages. The same gages were used to measure 

the lift force or the in-line force by rotating the gages 

90 degrees. It is easy to show that the calibration factor 

of the lift force (here in the vertical direction, up or 

down), expressed in terms of pounds or millimeter deflection, 

is Sina times the calibration factor for the in-line force. 

This proved to be true experimentally also and demonstrated 

independently the validity of the measurement technique. 
Raines [7] carried out experiments with the following 


cylinders: 


90 degrees, 5.905 inch, rough and smooth (nominal 


fon 


Size 6 inches) 
60 degrees, 6.005 inch, rough 
45 degrees, 5.923 inch, rough and smooth 


45 degrees, 4.42 inch, rough and smooth (nominal 
size 4.5 inches) 


45 degrees, 3.00 inch, rough and smooth (nominal 


size 3 inches) 


The vresent experiments were carried out with the 


following cylinders: 


60 degrees, 5.965 inch, smooth (nominal size 
6 inches) 


60 degrees, 4.450 inch, smooth and rough 
(nominal size 4.5 inches) 


60 degrees, 2.985 inch, smooth and rough 
(nominal size 3 inches) 


Relative sand roughness of k/D = 1/100 was used for all 


rough cylinders. 


ae 





In addition, the foregoing data were complemented with 
data reported in [10]. These data, taken using 3 inch and 
4.5 inch nominal diameter cylinders at 90 degree yaw angle, 


were used in the analysis and comparison of the entire data. 
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III. DATA ACQUISITION AND PROCESSING 


Data acquisition was accomplished using an HP3052A 
automatic data acquisition system. The system was configured 
as described in [9]. 

The data acquisition system consisted of an HP3495A 
20 channel scanner and an HP3437A system voltmeter, both 
controlled by an HP9845B desk top computer. The in-line 
or transverse signal, Originating as a voltage from the 
strain gages, was amplified and sent to the strip chart 
recorder. This amplified signal was also sent to the scan- 
ner aS one channel of data. The flow amplitude signal 
Originated from a differential pressure transducer, and was 
Similarly amplified and sent to both the strip chart recorder 
and the scanner as another channel of data. Strip chart data 
waS maintained as a visual record of the flow amplitude and 
force. The scanner switched between the two channels at 
designated intervals when triggered by the voltmeter. The 
voltmeter read each one of these voltage values and trans- 
ferred them to the computer for immediate calculation of 
the governing parameters. 

At least six cycles of in-line force data, recorded 
approximately 70 times per cycle, were acquired with the 
System. The data were then averaged to give one cycle from 


which the governing parameters were calculated. Ten cycles 
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of transverse force data, recorded at the Same rate, were 
acquired with the system. The transverse data was not 
averaged prior to calculating the governing parameters. 

The use of the HP3052A system proved highly beneficial 
for a number of reasons. The calculated parameters were 
available on a real time basis, allowing for detailed study 
of anomalies in the force patterns. Taking data over a 
large number of cycles removes more of the randomness of 
the flow forces. This was particularly valuable for 


analyzing the transverse force data. 
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IV. GOVERNING PARAMETERS 


Data reduction for the in-line forces is based on 
Morison's equation using the normal component of velocity. 
The derivation of the governing parameters is described in 
detail in [3 and 9]. Here only a brief description of the 
parameters will be given. 


For a Sinusoidally varying flow, Eq. (1) mav be 


written: 
2F ss 
—, = = c_ sine - C,!cose|cose (2) 
2 K ~“m d 
oDU 
ma 
where the flow is represented by U = -U_cosé Wie oe he 


The Fourier averages of om and C., are given as [1]: 


d 


ZT 
F sine 
ae ao de (3) 
m 3 D 4 
q 0 p m 
and 
3 em EF cose 
m 
== f de (4) 
da 4 2 
0 oD Un 


where Fi, represents the measured force. 
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As discussed in [9], it would be beneficial to remove 
the cycle-to-cycle variations in the force. This was accom- 


plished by using the normalized RMS value of the actual force 


ee.) - Two new coefficients were defined from this: 
ea — 
Car ~ Sone i Car ~ eae (3) 
and 
* 
Car ~ Se orane , Cnr = Sr corer (6) 
Using Eq. (2), the rms value of the calculated force 


can be shown to be: 


(7) 





Transverse force data were represented by the normalized 


rms value of the measured lift force (C ). Previously, 


iLrms 
some of the lift data were analyzed [7] using a Fourier 
transform to obtain frequency information. In the present 
study, this effort was discontinued because the frequency 
information so far obtained was more than sufficient to 


draw general conclusions regarding the frequency content of 


mae jift force. 
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Vo EXPERIMENTAL RESULTS 


fee PRESENTATION OF THE RESULTS 


G 


The force-transfer coefficients (Ce Car Caen ae 


ar’ 


Cc EE? are presented in graphical form in terms of the 


iPem 
governing parameters K, 8, a, and k/D in Figs. 2 through 7. 
The new as well as the previously-obtained data [7 and 


10] have been cross-plotted in Figs. 8 through 13 in order 


to exhibit the effects of the governing parameters. 


Eee DLOCUSSION OF THE RESULTS 
1. Inertia and Drag Coefficients 

The first indication of the fact that the independence 
princivle is not applicable to yawed cylinders in oscillating 
flow is seen in Fig. 8a where e 1s plotted as a function of 
K for 8 = 4040, k/D = 0.00, and a = 90, 60, and 45 degrees. 
Evidently, had the independence principle been valid, the 
three sets of data corresponding to three yaw angles would 
have collapsed into a single set, thus rendering the om 
versus K relationship independent of the yaw angle (at least 
Within the range of the test parameters and experimental errors). 
This conclusion is more universal than it may appear at first 
Sight on the basis of Fig. 8a. In fact, Figs. 8b through 8f, 
presenting Cx versus K for both smocth and rough cylinders of 
6, 4.5, and 3 inch nominal diameters, show that C. varies with 


the yaw angle as it did in Fig. 8a. Thus, one may conclude 
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that the independence principle is not valid at least for the 
inertia coefficient. The recognition of this fact is impor- 
tant for two reasons. Firstly, the inertial force is only 
part of the total force. Secondly, the significance of the 
inertial force depends on K. Thus, the failure of the 
inertia coefficient to obey the independence principle may 
not by itself be very important if the inertial force is 
small relative to the total force for the particular case 
under consideration. Consequently, it 1S important to 
examine the variation of the drag coefficient and its 
relation to the total force and the independence principle. 


Figure 9b shows the variation of C, as a function of 


a 
K for 8 = 2345, k/D = 0.00, and a = 90, 60, and 45 degrees. 


Clearly, the variation of C., with K and the yaw angle is far 


d 
more complex than that of Ca Figures 9a through 9c, 


depicting C., versus K for the three smooth cylinders, show 


d 
that the independence principle does not hold true except at 
relatively large values of K. The implication of this result 
is that when the vortex motion about the body becomes ex- 
tremely turbulent and incoherent, as it would at high K 
values, the overall effect of yaw becomes less significant. 
It is also in this region that the force is drag dominated 
and the drag force may be calculated in terms of the force 
acting on the projected area. Though the inertia coefficient 


exhibits strong variation with yaw even at high K values, 


this variation is of academic interest only since the 
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Bemteribution of the inertial force to the total force is 
merener inSilgnificant. 

Roughness plays an important and, oddly enough, a 
unifying role on the drag coefficient. This is somewhat 
anticipated on the basis of the fact that roughness precip- 
itates earlier tranSition and greater incoherence in the 
flow about a yawed cylinder, thereby creating conditions 
more favorable to the independence principle. Figures 9d 
through 9f show that the independence principle for rough 
cylinders is almost valid with the exception of the drag- 
inertia dominated regime (8<K<20). 

It is clear from the foregoing discusSion of the 
inertia and drag coefficients that the normal force acting 
on a yawed cylinder is Significantly underestimated through 
the use of Morison's equation, independence principle, and 
the drag and inertia coefficients appropriate to the normal 
cylinder. Only a narrow range of K values (from about 9 to 
17), in the drag-inertia dominated regime, that the drag 
coefficient for a normal cylinder exceeds that for a yawed 
cylinder. However, the large increase in the inertia 
coefficient of the yawed cylinder in the same region of K 
values more than offsets the relatively small decrease in 
Mes drag coefficient. 

The absence of the inertia crisis, the large increase 
in the inertia coefficient, the relative decrease in the drag 


and lift coefficients, all in the drag-inertia dominated 
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regime, of yawed cylinders are a consequence of the decrease 
of the spanwise coherence of vortices. In general, the 
smaller the spanwise coherence, the larger the Ca Eh eines 
region. It appears that the independence principle or the 
cosine law, as it 1S sometimes called, is a gross simplifica- 
tion of the behavior of flow in the near wake. One may, 
therefore, conclude that the Fourler-averaged drag and inertia 
coefficients, based on Morison's equation, are unique for 
each angle of yaw, Reynolds number, Keulegan-Carpenter number, 
and the relative roughness. 

2. RMS Values of the Inertia and Drag Coefficients 

It has already been shown [7] that the normalized rms 

values of the inertia or drag coefficients remove signifi- 
cantly the scatter from the data and establish clearer 
relationships between these coefficients and the governing 
parameters. The same rms coefficients are used here to 
examine the effect of yaw and experimental error. Figures 
lla through 1l1f show C= = c_/C asa Lupetwon ofek. 


arms 


Seearly, eae is independent of Reynolds number (or 8), yaw 


angle, and the roughness for K smaller than about 8. For 
smooth cylinders, on is nearly independent of yaw for K 
larger than about 20. Additionally, for rough cylinders Car 
is least dependent on yaw for relatively low values of Re. 


Whe drag coefficient C exhibits similar character- 


Oule 


istics (see Figs. l2a through 12f). The scatter in Car is 


considerably less than that in C., primarily because of the 


d 
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fact that the variations in Ca aS a consequence of the random 
disturbances imposed on the measured force are smoothened out 
by the variations in the corresponding ae values. This 


fact has been substantiated through numerous examples by 


examining the relative variations of Car Cat ana: € As 


arms ~ 


K increases C reaches its theoretical value of 1.63, as 


aie 
seen from Eq. (7). Finally, roughness tends to decrease the 
dependency of Car on Re (see Figs. 12d through 12f). 

3. Transverse (Lift) Force 

The lift force data are presented in terms of the 

normalized rms values of the measured lift force in Figs. l3a 
through 1l3c for the three smooth cylinders. These figures 
show that Ces 1s largest for the 90-degree cylinders, as 


one would expect, and decreases rapidly with yaw. However, 


meme VYarlation of C with the yaw angle (60 and 45 degrees) 


irms 
1S not as large as that with the Reynolds number. The lift 
coefficient increases rapidly as K increases from about 6 


memi2. Then, C again decreases rapidly due to the 


lrms 
increased lack of coherence of the vortices. Similar 
observations may be made for rough cylinders (see Figs. 13d 
through 13f). It suffices to note that the yaw has an 
immediate effect onthe coherence of vortices and on the lift 
coefficient, particularly in the drag-inertia dominated 
regime. For large values of K or Re, the irregularity of 


the vortex shedding even for a 90-degree cylinder obscures 


the yaw effect. It is then rather difficult to distinguish 
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between the lift coefficients of cylinders at various yaw 
angles. Inthe ocean environment, one is more likely to find 

a confused vortex shedding with very small coherence length. 
Thus, One may safely state that the inclined members subjected 
to the wave motion in the ocean environment will experience 


considerably less lift force than those normal to the waves. 


Bn3 





Vals DENS EPpy et YY ANALYSES 


Pee PHASE SHIFT EFFECTS 

It has been known Since the inception of Morison's 
equation in 1950 that the phase change between the maximum 
force and the maximum velocity may lead to varying degrees 
of error, depending on the governing parameters. This phase 
shift may arise due to electronic or numerical filtering of 
the force and velocity data or due to variations in the 
coherence of the vortices along a given separation line. 
Even though it 1S not yet possible to attribute a given 
fraction of the estimated error to a Specific cause, it is 
most desirable to examine in detail the isolated effects of 
such things as an imposed phase shift and/or a random 
disturbance of a given amplitude and frequency. Eventually, 
such information will show the region of sensitivity of the 
force coefficients and Morison's equation to the disturbances 
imposed by the environmental conditions and will lead to a 
more rational interpretation of the causes of large scatter 
observed in the results of ocean tests. 

The effects of phase shift onthe force-transfer 
coefficients have been examined as follows. For a given 
K, 8, k/D, and yaw angle the drag and inertia coefficients 
Meee ocen Calculated from Eqs. (3) and (4). These 


coefficients correspond to cases where there was no phase 
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shift as ascertained by the manner in wnich the experiments 
were carried out. Then the force was recalculated through 
the use of Morison's equation and the above coefficients. 
The calculated force was then given a phase shift and a new 
set of force coefficients were calculated through the use of 
Eqs. (3) and (4). The method was repeated for representa- 
tive values of K and reasonable phase angles. 

Figures 14a through 14e show the results for five 
coefficients. Clearly, the phase angle has significant 
effect on all of these coefficients with the excepotion of 
ens. Concentrating attention on Cs and Ca it is noted 
that oN decreases and Ca increases with increasing phase 
angle (where the force curve is shifted forward in time). 
The relative rate of decrease or increase of these two 
coefficients certainly depends onthe particular values of 
the governing parameters. For example, for K = 6.66 (Fig. 
14a) om decreases rather gradually whereas Ca increases 
rapidly with increasing phase angle. The reason for this 
is that for values of K less than about 8 the fluid force 
is inertia dominated and the inertia coefficient is not 
very sensitive to small disturbances imposed on it. In 
mathematical terms the foregoing simply means that the force 
curve is like a sine curve and that when it is shifted right 
or left by some amount and then multiplied with sing, as in 
Eq. (3), the integral could not be too sensitive to the 


phase difference. The drag coefficient, on the other hand, 
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is extremely sensitive to small phase shifts in the inertia 
dominated region simply because a small phase shift ina 
Weane=-like" curve multiplied with |[cos@|cos9, as in Eq. (4), 
produces large changes. 

In the drag-inertia dominated regime (K ranging from 


meout 8 to 20) both C. and Cn are affected, as seen in Figs. 


d 


l4c and 14d. In the drag dominated regime, however, Ca 1s 
hardly affected by phase shifts; whereas om shows a dramatic 
drop with increasing phase angle (Fig. 14e). The reason for 
this 1s easily explainable in a manner similar to that 
previously presented in connection with Fig. 14a. 

The overall conclusion to be drawn from the foregoing 
is that there is no flow regime where one or the other 
coefficients is not effected by the phase shift. In fact, 
if a constant phase shift of +5 degrees were to be randomly 
imposed on a series of data obtained with K values ranging 
from 5 to 100 one would have observed large scatter in Cy 
and small scatter in oe for K smaller than about 8; large 
Scatter in both C4 and Coa bourse Ecom apout o to 20> and 
large scatter in Ca and small scatter in C4 for K larger 
than about 20. These scatter patterns are invariably 
observed in the data obtained in the ocean environment. 
As noted earlier, not all of the phase shift is due to 
electronic filtering of the data. Some of the phase shift 
comes from the reduction of the coherence length of the 


vortices by the disturbances ever-present in the wave field. 
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Thus the purpose of the foregoing is not to suggest that 

all phase shifts could be avoided or eliminated but rather 

to quantify the effects of the phase shift in terms of the 
expected variation of the force coefficients in order to have 


a better grasp of the meaning of the observed scatter. 


B. RANDOM SQUARE WAVE EFFECTS 

The experiments conducted under controlled laboratory 
conditions are immune to practically all of the disturbances 
encountered in the ocean environment. Consequently, the 
vortices have larger coherence lengths and these give rise 
to maximum values of the semi-peak-to-peak force coefficient 
C. (spp) .* In the ocean environment, however, the variable 
wave direction, change of wave velocity with depth, effect 
of secondary waves with all sorts of harmonics, irregular 
surface roughness, and numerous other stochastic variables 
give rise to varying degrees of scatter in the force 
coefficients based on ocean data. It 1S important from the 
engineering point of view to examine the collective impact 
of these disturbances on the force coefficients. It 1s under- 
stood that very high frequency disturbances imposed on the 
measured force cannot significantly affect Ca and C.. On 


the other hand, disturbances which simulate either in part 


* 

This does not imply that perfect coherence leads to 
ieeimum values Of CG and Cy. In fact, a force trace 
resulting from a perfect coherence can be rendered 
incoherent by a random square wave. The new force 
trace yields a larger Ce and a smaller Ca: 
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Or in combination a "Sine like" or a “cosine like" wave are 
most effective in causing large changes in Ca and Ca ies 
follows from the very definition of the inertia and drag 
coefficients as in Egs. (3) and (4). Clearly one must 
examine the effect of randomly imposed disturbances of 

given frequency and phase shift on the variation of the 

force coefficients. Such an extensive analysis has been 
carried out in a manner described as follows. A frequency 
mat1o of E ug/ Ey was chosen where oe is the frequency of 

the square wave and f, is the frequency of the periodic 
oscillation of the flow in the tunnel. Then the square wave 
was asSigned a random amplitude and initial phase shift 
through the use of a built-in random number generator. The 
maximum of the disturbance was taken to be the absolute of 
C,(spp) (almost equal to the absolute maximum Of the measured 
force). This procedure varied the amplitude of each disturb- 
ance by multiplying a random number between +1 and -1 with 
Cp(spp). Calculations have been performed 1000 times for a 
given frequency ratio and force trace. Each calculation 
started with a randomly selected phase angle and followed 
with square wave fluctuations whose amplitude varied randomly 
as noted above and whose frequency was fixed by the frequency 
ratio. Figures 15a through 1l5e show sample plots of the 
Original force, the disturbed force, and the square wave 
disturbance. These figures are presented here for the 


purpose of illustration and not necessarily for the 
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@mantification of the boundaries of the sensitivity 
envelope. 

Figure 16a shows the results obtained in the manner 
described previously. The two symbols in this figure 
represent the relative changes in Cn and Ca: The computer 
program was directed to identify the particular cases in 
1000 runs which gave the largest changes in Ca and the 
corresponding changes in Ca: It is important to emphasize 
that the change in Ca in Fig. l6a is not the maximum 
possible change that could have taken place with a suitable 
Square wave disturbance. Thus, Fig. 16a shows that there 
are certain frequency ratios which cause a large change in 
ern with very little change in the corresponding Ca: As 
expected, disturbances with a frequency ratio of 1.0 yield 
the largest AC _/C, Since this nearly corresponds to a sine 
function. Surprisingly enough, high frequency disturbances 
have very little effect on the inertia coefficient, at least 
in the drag-inertia dominated regime (see Fig. l6a). fThus, 
one can filter out such high frequency disturbances without 
materially affecting the determination of the force-transfer 
coefficients. 

A similar calculation has been performed with the force 
trace corresponding to Fig. 16a to determine the maximum 
change in Cy and the corresponding change in Ca LOY 
representative frequency ratios. Figure 16b shows the 


results obtained in this manner. Several facts are 
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immediately apparent. C, 1s not as sensitive as Cc, to the 
randomly imposed disturbances in the drag-inertia regime 
even though AC3/Ca is guite appreciable. This is directly 
related to the stability of the integral representing Car 
i.e., |cos@|cos6 times asquare wave. 

Figures 17a and 17b show the effect of the disturbances 
on om and Cy in the inertia dominated regime and Figs. 18a 
ema 18b illustrate the similar effects in the drag-inertia 
dominated regime. It is clear from all of the figures 
pertaining to the sensitivity of the drag and inertia 
coefficients that it 1s not the dominance of the drag or 
inertia in a given regime that determines the sensitivity of 
eee emecoctiticient. To be sure, the said effect is present 
in Figs. l6a through 18b. However, the most significant 
elements of a disturbance are its frequency and relative 
magnitude. Disturbances imvosed on the force at a frequency 
ratio of about one or unity will cause major changes in 
both Cy and Cn 1rrespective of the flow regime. As to the 
amplitude of the disturbances, they may depend on a number 
of parameters such as K, Re, turbulence in ambient flow, 
interference of other bodies, etc. In the present analysis 
the maximum amplitude of the disturbances at any frequency 
were taken to be equal to C,,(Spp) of the measured force. 
Since this coefficient decreases with increasing K, it 1S 


easy to understand that at high K values AC Cm oe AC./C, 
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will result in values as large as those encountered at 
smaller K values (see Figs. 16a through 18b). 

The significance of these results will be illustrated 
with a numerical example. Assume K = 14.08, Beta = 4040 
fis in Fig. l6a), for which Cn = 9-61, Cy = 1.63 and C(spp) = 


1.66. For a frequency ratio of Be) ae = 1.0, one has 


fe /C_. 


Assuming a disturbance amplitude equal to 10 percent of 


2. 7] ame AC 3/C 4 =e? Letsomerrcs. lGa and lob. 


C,(spp), one has: 


Sepeaturped §_ in = . 
d Sa es = (LEAC./C, * 3C,(Spp) ) * C4 (original 
value) 
1.78 
= a * = 


and 


C (perturbed _ A a r: hone 
eA Lue) Se. Cae 3C,(Spp) ) i 


at ass ° 
Cee OL 0.6) = 0.45 


The foregoing examples show that a random disturbance of 
frequency ratio 1.0, imposed on a particular force, may 
have resulted in relatively large cycle to cycle variations 


of the force coefficients. Here Ca Varies ErOM a minimum 
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value of 1.48 to a maximum value of 1.78 and Ch varies from 

a minimum value of 0.45 to a maximum value of 0.77. These 
are relatively large changes and show clearly the importance 
of the following two factors. First, it 1s not possible 

to obtain scatter-free data in the ocean environment even 

for a Single K and Re because of the ever-present random 
disturbances emanating from various sources. Secondly, a 
structure cannot be designed to withstand the maximum forces 
predicted on the basis of idealized laboratory experiments 
because the maximum forces may occur only in a small fraction 
of the lifespan of the structure. Clearly, beyond this point 
one has to make subjective decisions regarding the probability 
of the occurrence Of maximum force. Figures Similar to l6a 
through 18b enable one to predict the force coefficients 
which would have prevailed under less favorable conditions 
(reduced coherence, directional waves, body interference, 


etc.). 
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Vr econe uu snONS 


The forces acting on yawed smooth and rough circular 
cylinders in sinusoidally oscillating uniform flow have been 
extensively investigated, and the following conclusions have 
been reached: 

1. The independence principle does not apply over signif- 
icant portions of the range of K values investigated. The 
Mimag and inertia coefficients differ significantly from those 
deduced from the 90-degree case for all 45-degree and 60- 
degree smooth and rough cylinder cases. 

2. The effect of roughness is to bring the yawed and 
unyawed cylinder force-coefficients closer, except in the 
drag-inertia dominated regime. 

3. Yaw has a very strong effect on the coherence of 
vortices and hence on the transverse force coefficient, 
particularly in the drag-inertia dominated regime. The lift 
coefficient is largest for the 90-degree cylinder, as might 
have been expected, and decreases sharply with yaw. 

4, An extenSive sensitivity analysis has shown that the 
force-transfer coefficients are affected by the phase shift 
between force and velocity to warious amounts dependent on 
the flow regime and the magnitude of the phase shift. 

5. Finally, the results have shown that the force 


coefficients and the validity of the Morison equation are 
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sensitive to disturbances imposed on the’ in-line force, 
particularly in the frequency range where the disturbance 
frequency matches the wave frequency. The stability envelope 
has been delineated through extensive calculations to enable 


one to interpret the scatter of the data obtained in the 


Beean environment. 
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